Creating a supersolid in one-dimensional Bose mixtures 
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We identify a one-dimensional supersolid phase in a binary mixture of near-hardcore bosons 
with weak, local interspecies repulsion. We find realistic conditions under which such a phase, 
defined here as the coexistence of quasi-superfluidity and quasi-charge density wave order, can be 
produced and observed in finite ultra-cold atom systems in a harmonic trap. Our analysis is based 
on Luttinger liquid theory supported with numerical calculations using the time-evolving block 
decimation method. Clear experimental signatures of these two orders can be found, respectively, in 
time-of-flight interference patterns, and the structure factor S(k) derived from density correlations. 
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The supersolid (SS) phase, defined as a many-body 
state that simultaneously shows superfluid (SF) and 
charge density wave (CDW), i.e. crystalline, order, has 
been an intriguing notion since its first proposal [l[, due 
to its seemingly paradoxical nature. Numerous studies 
of SS phases [2j, |3[ have recently been reported, moti- 
vated by fundamental theoretical interest in a system 
that exhibits competing orders, and by recent experimen- 
tal reports of observations of supersolidity of 4 He in vycor 
glass 0, [B[ . This 4 He system exemplifies the complexity 
of studying strongly correlated systems in a solid-state 
context: it combines strong disorder due to the porous 
medium, strong interactions between atoms and of atoms 
with surfaces. Under such circumstances it is difficult to 
demonstrate the existence of a SS phase, which involves 
a subtle competition of fluctuations. 

In this paper we show that supersolidity can be stud- 
ied with clarity in another physical system: ultra-cold 
atoms in optical lattices. Since the demonstration of the 
SF-Mott insulator transition in 3D , the technology of 
cooling and trapping atoms has supported studies of nu- 
merous quantum many-body phenomena, such as BEC- 
BCS crossover Q, noise correlations [§], the Berezinsky- 
Kosterlitz-Thouless transition the Tonks-Girardeau 



gas 



lfl llj . transport and collisional properties of one- 
dimensional (ID) gases O, and the Mott transition in 



ID [lH and 2D 14 |. Appealing features of this technol- 
ogy, from the perspective of many-body theory, are that 
it creates well defined and tunable systems, and that the 
set of measurable quantities differs from those in solid- 
state systems. Thus, ultra-cold atom systems can give 
interesting and unusual insights into many-body states. 

The objective of this paper is to propose a realistic 
setup of how to create and detect a supersolid with cur- 
rent technology. Specifically, a binary mixture of near- 
hardcore bosons with weakly repulsive inter-species con- 
tact interactions in a ID potential displays both CDW 
and SF quasi- long range order (QLRO). Such mixtures 
have an inherent tendency to undergo phase separation, 
which can be avoided if the inter-species interactions are 



sufficiently weak. 

We study the SS phase with analytical and numerical 
techniques, and present a concrete proposal for its real- 
ization in current experimental systems. First we use a 
Luttinger liquid (LL) approach to derive the phase dia- 
gram of the homogeneous, infinite system, with a renor- 
malization group (RG) calculation. We then address the 
question of realizing such a phase under realistic condi- 
tions, for a finite system of ~ 10 2 lattice sites in a har- 
monic trap. Using a number-conserving time-evolving 
block decimation (TEBD) method [HI we numerically 
determine, with a well-controlled error, the ground state 
of the system from which we extract various correlation 
functions. We first identify the SS phase, through signa- 
tures in the pair and anti-pair correlations. Other cor- 
relations contain information that is accessible to direct 
experimental observation, and we discuss possible exper- 
imental signatures of the SS phase, i.e. the coexistence 
of SF and CDW order; the SF order is manifest in the 
single-particle correlation function, which can be deter- 
mined from time-of-flight (TOF) interference patterns; 
the CDW order is seen in density-density correlations, 
which is reflected in a measurable structure factor. 

We consider a mixture of two species of bosonic atoms 
with short-range interparticle interactions, confined in a 
ID optical lattice and an additional harmonic potential. 
Contemporary experimental realizations of such systems 
are usually well approximated by a Hubbard model: 

H = -t ]T bl A b aJ + — ^2n ati (n a ,i-l) 

<ij> ,a i , a 

+ Uu ^ n l,i n 2,t + X! ^i 2 ( n l,j + n 2j)- (i) 



Here t is the hopping energy; b a ,i is a boson field operator, 
with a = 1,2 a species index and i a lattice site index; 
U (U12) is the intra- (inter-)species on-site interaction 
energy; n a> i = b' ai b a ^; and f2 represents the strength of 
the harmonic trap, which is centered on the site j = 0. 
We now derive the phase diagram of this system from 
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FIG. 1: Phase diagram of a Bose mixture, from LL theory, as 
a function of the inter-species interaction U12 (in units of v, 
see text), and the LL parameter K of the uncoupled system 
(Eq. ©); on the left the full diagram, on the right (cor- 
responding to the shaded area) the near-hardcore, repulsive 
regime with the supersolid (SS) phase. For attractive inter- 
actions, in the paired regime, paired SF (PSF) is dominant, 
with CDW QLRO being subdominant in parts of it. Out- 
side of that regime, SF is the dominant quasi-order, but for 
the repulsive, near-hardcore regime we have CDW QLRO as 
well, which constitutes a supersolid phase. For large repul- 
sive interactions the system phase separates (PS), for large 
attractive ones it collapses (CL). 



LL theory. We consider two ID bosonic SFs, with den- 
sities equal to each other, but incommensurate to the 
optical lattice. In particular wc exclude half- and unit- 
filling, which would destroy the SF order, by choosing the 
density and the global trap in such a way that even at the 
trap center the density of each species stays below 0.5. 
The essential function of the optical lattice is to provide 
a sufficiently large ratio oiU/t, as in [IcJ ■ We emphasize 
that the supersolid phase also exists in the absence of 
a lattice[3]; however, since present quasi- homogeneous 
ultracold atomic systems are generically closer to weak 
coupling, the presence of a lattice is advantageous to ex- 
perimental realization of supersolidity. 

The basic concept of LL theory is to express the 
bosonic operators b a ,i through a bosonization identity, 
such as Haldane's construction [T7L fl8l] : 

b h2 (x) = [n + n 1 , 2 (a;)] 1/2 ^e 2miei ' 2(a;) e l * 1 ' 2(3;) ,(2) 



where we switched to a continuum model, b a ^ —* b a (x). 
n is the average density of the two species, Hi^(x) are 
the low-k parts (i.e. k <C 1/n) of the density fluctuations; 
the fields 81,2(2;) are given by 0^2(2;) = nnx + 0i^(x), 
with #1,2(2:) = 7r J x dyLIi^y). $1,2 are the phase fields, 
the conjugate fields of the density fluctuations LTi^x). 

In terms of these fields, the action of the two coupled 
bosonic SFs is given by [3, [2(| : 

J 3 = 1,2 




FIG. 2: The single-atom density of species 1 (red) and 2 (blue) 
of a system of 90 sites, with 17 atoms of each type, with t/U = 
0.005, U12/U = 0.04, and a trap parameter Q = 10~ 5 (7. 
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cos(26»i - 20 2 ) (3) 



The two SFs are characterized by a LL parameter K and 
a velocity v, which is contained in r = (vt,x). The LL 
parameter K is a measure of the intra-species interaction; 
in the near-hardcore regime we use [19j : 
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(4) 



Similarly, v can be related to the parameters of the un- 
derlying Hubbard model by v ~ vf(1 — 8tn cos nn/U), 
where vf is the 'Fermi velocity' of an identical system of 
fermions, vp = 2tsin7rn. The density-density interaction 
between the two SFs creates both the term containing 
V#iV#2) as well as the backscattering term 0, [S^, con- 
taining cos(2$i — 2#2)j which describes short-range inter- 
species repulsion. We change variables to the symmetric 
and antisymmetric combinations 4> s / a = -^=(<f>i±(f>2) and 
9 s / a = -^(0i ±^2), and diagonalize the quadratic part of 
the action which gives the following parameters for the 
two sectors: 



K s/a = (1/K 2 ±U 12 /vttK) 



-1/2 



(•5) 



which to lowest order gives K s / a K^fU^K 2 /2irv. The 
effective velocities are v 8 / a = ^v 2 ± U 12 Kv/n. Phase 
separation (collapse) is reached when v a r s ) becomes imag- 
inary. The anti-symmetric sector contains the nonlinear 
backscattering term. To study its effect, we use an RG 
approach; the flow equations for which are given by[20j: 



dl 



(2 - 2K a )g a 



dK a 
dl 



a 2 

0_2 a 



(6) 



This set of flow equations has two qualitatively different 
fixed points: Either g a diverges, driving a pairing tran- 
sition, which in turn renormalizes K a to zero, or g a is 
renormalized to zero. In the latter case the Gaussian 
fixed point is restored with a finite effective value K*. 
Therefore the correlation functions are again algebraic, 
containing this effective parameter. It is this second sce- 
nario that we are interested in, not the actual phase tran- 
sition itself. 
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FIG. 3: LL parameters K a and K a , as a function of U12, for 
the same parameters as in Fig. [2] from numerical fitting. 



To determine the phase diagram we consider the corre- 
lation functions of these order parameters: single-particle 
SF, described by Osf = b a , CDW order, corresponding 
to the 2fcF-component of the density operator Ocdw = 
n a , and paired SF, Opsf = &i^2, which appears on the 
attractive side. The form of these correlation functions 
is (O(x)O(O)) ~ |x| Q_2 , except for the single-particle SF 
in the paired regime, where it decays exponentially. An 
order parameter O(x) has QLRO, if its correlation func- 
tion is algebraic, and a > 0. This implies that the cor- 
responding susceptibility is divergent, indicating an in- 
stability towards ordering [2(j. The scaling exponent of 
Osf is «sf = 2 — l/AK s — l/4K a , the one of Ocdw is 
acDW = 2 — K s ~ K a , and PSF has apsF = 2 — 1/K S . 
We also consider the anti-pair operator b\ &2 , which has a 
scaling exponent of 2 — 1/K a . We use the latter two cor- 
relation functions in the numerical fitting procedure. In 
Fig. [T]we see the resulting phase diagram. For attractive 
interactions we see the formation of a paired phase, in 
which two regimes of quasi-order are found: in the entire 
paired regime, PSF is the dominant QLRO, whereas for 
part of that regime we find CDW as a subdominant or- 
der. The latter can be considered a SS of pairs, whereas 
single-particle SF is destroyed. On the repulsive side we 
find the SS phase that we look for in this paper. Using 
the flow invariant g 2 — An 2 (K a — l) 2 , and Eq. (j4|), we 
determine the nearly linear SS phase boundary for small 
repulsive interactions to be: Unjv > 32t sunrn/U. 

Having given the phase diagram of the infinite, homo- 
geneous system, we now address the question of how su- 
persolidity can be found in actual cold atom systems. For 
this, we use a TEBD method 15[ to obtain the ground 
state of the system. We choose a small value for t/U, 
to reach the near-hardcore regime, and a positive U12 of 
the order of t, to avoid phase separation. We choose the 
atom number and the global trap parameter, such that 
the density is smaller than 0.5 throughout the system. 
In Fig. [2] we show the densities of the two species for 
the case t/U = 0.005, U 12 /U = 0.04, fl = 10~ 5 [/, and a 
particle number of 17 of each atom species, on a lattice of 
90 sites, as an illustration of the ground state. One can 
clearly see the density modulation of each species, whose 
wavelength is determined by the density, not the lattice. 




FIG. 4: Interference pattern of the atomic mixture, with the 
parameters of Fig. [2] when interpreted as 87 Rb, released from 
an optical lattice with lattice constant a — 400nm, localized 
states with a length scale d = 60nm, after an expansion time 
t = 30ms. This could be realized with an optical lattice po- 
tential of 18Er in longitudinal, and 30_Er in transversal direc- 
tion, En being the recoil energy. In a) we see the ID density 
n of the full pattern (in units mm -1 ), as a function of the spa- 
tial coordinate x (in mm). In b) we fit the central interference 
peak with a power-law. The good numerical fit indicates the 
presence of single-particle quasi-superfluidity. 



A central question when studying 'phases' in finite-size, 
non-homogeneous systems, is whether a given state can 
be reasonably related to a phase of the associated system 
in the thermodynamic limit. We address this question by 
numerically fitting the correlation functions of the pair 
and the anti-pair operator in the bulk of the system with 
power-law functions. We find a very good fit, and we de- 
pict the scaling exponents extracted from the numerical 
data in Fig. [31 as a function of the interaction C/ 12 . 

Having established that the state of the system is in- 
deed a supersolid, we now turn to the crucial question 
of how this phase can be detected in experiment. We 
propose two measurements that would address this ques- 
tion: 1) a TOF interference measurement to determine 
the single-particle correlation function, which is the defin- 
ing quantity of SF QLRO, and 2) a measurement of the 
structure factor [21( to determine the density correlation 
function, which is the defining quantity of CDW QLRO. 

TOF measurement. We assume that when the op- 
tical lattice is turned off, the atoms expand freely, 
that is b a (x 1 t) = J2j w ( x ~ r j)ba,j, where w(x,t) — 

y / d/\/2¥A(t) 2 exp(-a; 2 /4A(i) 2 ), with A(t) 2 = d 2 + 
ith/2m, a the lattice constant, d the width of the ini- 
tial state, assumed Gaussian, t the expansion time, and 
m the atomic mass. We calculate the density n(x) — 
(b^ (x,t)b(x,t)) , which contains the single particle corre- 
lation function of the original system, which we show in 
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FIG. 5: Structure factor of the atomic mixture, with the 
same parameters as in Fig. 2] in 1/fim, as a function of k 
in units of 2n/a. a) corresponds to very weak interactions 
Un/U — 0.0025, and shows no particular structure beyond a 
SF signature, b) corresponds to U12/U = 0.04. We see addi- 
tional peaks, that indicate CDW QLRO at a momentum that 
is consistent with the density at the center of the trap. 

Fig. 2J Fig. Q] a) shows the full interference pattern, and 
2] b) the central peak with a fit with a power-law func- 
tion c(x 2 + a 2 ) / 2 ' 1 , shown on a log-log scale. The good 
agreement with power-law scaling indicates the presence 
of SF QLRO. 

Structure factor. We assume, as an example, an in- 
situ measurement of the same system described before. 
The real-space structure factor is related to the density- 
density correlation function of the lattice system by: 

S (k) « J^l e - ik ^({n xl n X2 )-{n xl ){n X2 ))(7) 

with X12 = x\ — X2, and the real-space correlation func- 
tion defined as 

(n Xl n X2 ) = 2J \w(xi - r il )\ 2 \w(x 2 - r i2 )\ 2 (n^n^) . 

(n Xl )(n X2 ) is similarly defined. Fig. [5] shows S(k) for 
U12/U = 0.04, and a nearly non-interacting example, 
U12/U = 0.0025. The envelope of the function is given by 
the inverse width of the Wannier state; the periodic shape 
comes about because we map a lattice quantity onto real 
space. We clearly see the onset of a peaked structure 
at momenta that are consistent with the density at the 
trap center. These peaks arc related to algebraic cusps 
in the static structure factor, when the CDW regime is 
reached, S(2kp + q) ~ |g| Q co»'-i a measurement of 
the dynamic structure factor would lead to an even more 
striking result, since these cusps would translate into al- 
gebraic divergencies. We also note that a measurement 
of the lattice density in-situ[22j could depict a profile as 
in Fig. [21 and that repeated measurements and a cor- 
relation analysis as in |8j, but without expansion, could 
give the real-space density correlation function directly. 

In conclusion, we have presented a proposal of how a 
supersolid phase of binary bosonic mixtures in ID can 



be created and probed in present ultracold atom exper- 
iments: in particular, with only local interatomic in- 
teractions. Using LL theory, we identified the generic 
phase diagram of this system for incommensurate filling; 
with TEBD simulations we found a concrete example of 
a finite, realistic system, including a global trap, which 
shows both SF and CDW QLRO. Two well-established 
measurement techniques, the TOF signal and the struc- 
ture factor, provide clear experimental signatures of the 
two orders present in this remarkable state of matter. 
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